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Abstract 

We approximate the Bolker-Pacala model of population dynamics 
with the logistic Markov chain and analyze the latter. We find the 
asymptotics of the degenerated hypergeometric function and use these 
to prove a local CLT and large deviations result. We also state global 
limit theorems and obtain asymptotics for the first passage time to 
the boundary of a large interval. 



1 Introduction 



The central problem in population dynamics is the construction of "realistic" 
models with non-trivial ergodic limits, that is, to construct birth-and-death 

*Thc first author was supported in part by NSF grant DMS-0739164. 
Key words Bolker-Pacala model; logistic Markov chain; limit theorems 
AMS Subject Classification 60J80,92D25 



1 



Markov processes that tend in law to a stationary state. These models must 
satisfy, in addition, several other requirements. First, they must possess 
stability with respect to random noise, which always exists in the environ- 
ment. Second, they should contain stochastic dynamics (we have developed, 
for example, models involving migration and, in some cases, immigration 
[6]). Third, because all bio-populations demonstrate an intermittency effect 
or "patches," these models should exhibit deviations from local Poissonian 
(uniform) statistics. 

The classical Galton- Watson process (see [9, 11]), for example, cannot 
satisfy these requirements, as it has two significant defects: there are no 
spatial dynamics - indeed, no spatial distribution of particles - and there is 
no interaction between particles. 

The first defect can be easily remedied if we have the particles undergo 
independent random walks. Recent work [6, 14, 16] contains the important 
result that if the underlying random walk is transient, then, under mild 
technical conditions, the point field n(t,x) has a unique stationary ergodic 
limit n{oo,x). Unfortunately, situations with other requirements are much 
more complex. 

The fact that the population is stable only in the critical condition that 
/3 = /i is less than satisfactory, given that, in the short run at least, there is 
no obvious reason why this condition should hold. Another troublesome fact 
arises from the work of Kondratiev, Kutoviy, and Molchanov [14]. They show 
that for a population in a time-independent, stationary-in-space random en- 
vironment, there is no non-trivial limiting distribution. Roughly speaking, 
then, the population is stable only if the parameters specify the critical con- 
dition at almost every point in the space. Clearly, it would be preferable to 
have an alternative model that more robustly yielded a stable distribution- 
which, after all, is found frequently in nature. 

For this we turn to the Bolker-Pacala model [1, 2]. The Bolker-Pacala 
model, well-known in the theory of population dynamics, is a stochastic 
spatial model that incorporates both spatial dynamics and competition. It 
is of especial concern, however, that we have no good qualitative information 
about the B-P (or logistic) process. Consequently, in this paper we begin 
to study the B-P model. Our approach is to consider a discrete version on 
the lattice Z*^ and then to use a random walk representation in mean field 
approximation . 

The paper is organized as follows. Section 2 offers a detailed descrip- 
tion of the Bolker-Pacala process. The next section describes our mean-field 
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model and presents the logistic Markov chain, which we use to analyze the 
process. Section 4 contains our analysis of the logistic Markov chain. These 
include, after preliminary supporting sections, limit theorems for the invari- 
ant distribution, global limit theorems, and analysis of the first passage time 
to the boundary points of an interval. 

Our principal results are contained in section 4. We begin by deriving 
asymptotic properties of the degenerated (confluent) hypergeometric func- 
tion, which we use to prove a local Central Limit Theorem. Subsequently, 
we obtain a large deviations result for the invariant distribution for the logis- 
tic Markov chain. Namely, if we make the weight of the quadratic term in the 
B-P model proportional to then the invariant probability for large de- 
viations from the equilibrium state is dominated by times a constant 
times the exponential of —L (Theorem 4.4). Based on theorems of Kurtz 
[17, 18], we state global limit theorems — a functional LLN and CLT — for the 
logistic Markov chain. Lastly, we obtain asymptotics for the first passage 
time from the equilibrium point of the logistic Markov chain. Specifically, we 
show that the first passage time from the equilibrium point to the bounds of 
a symmetric interval, the left bound of which is close to 0, is dominated by 
L~^/^ times a constant times the exponential of L (equation 4.23). 

2 Preliminaries: Description of the Process 

In this section, we introduce the general Bolker-Pacala model, which can 
be formulated as follows. At time t = 0, we have an initial homogeneous 
population in M"^, (i > 1, that is, a locally finite random point process 

no(r) = ^(particles in F at time t = 0). 

The simplest option is for no(r) to be a Poissonian point field with intensity 
Po > 0, i.e., 

P{n,{V) = k} = exp(-p|r|)^^lU^, A; = 0, 1, 2, . . . 

The following rules dictate the evolution of the field: 

i) Each particle, independent of the others, during time interval {t,t + dt) 
can produce a new particle (offspring or seed) with probability A'^dt, 
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A'^ > 0. The initial particle remains at its initial position x but the 
offspring jumps to x + z + dz with probability 



(We stress that this differs from the classical branching process, in 
which the "parental" particle and its offspring commence independent 
motion from the same point.) We assume, of course, that all offspring 
evolve independently according to the same rules. 

ii) Each particle at point x during the time interval (t, t + dt) dies with 
probability /irft, where yU. is the mortality rate. 

iii) Most important is the competition factor. If two particles are located at 
the points x,y & W^, then each of them dies with probability a~{x—y)dt 
during the time interval {t,t + dt). This requires, of course, that a~{-) 
be integrable; set 



Here we have interacting particles, in contrast to the usual branching process. 
One can expect physically that for arbitrary non-trivial competition (a~ G 
(^(R*^), A~ > 0), there exists a limiting distribution of the particles. This 
is expected to be true because, heuristically, if the population is small at 
each site, the competition will also be very small, but b > fi will allow the 
population to grow. On the other hand, if the population is very large, 
competition is very large, and death is proportional to n?{t,x), due to the 
competition term, while birth is proportional to n{t,x), thus the population 
should decrease. At present, this fact has been proven only under strong 
restrictions on a"*" and a~ [7, 8]. 



In this section, we introduce the mean field approximation to the general 
Bolker-Pacala process. Let C Z'^ be a box with \Ql\ = L, L being a large 
parameter, and suppose that no particles exist outside of Ql. We let 




A 




3 Mean Field Approximation 



a~^{x) = — on Ql, a (x) 



T? on Ql 
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for some rates k, 7 > 0. Thus, after a jump each particle is distributed uni- 
formly on Ql. We let h and /i be the birth and mortality rates, respectively. 
Let 

NL{t) = nit,x) 

be the total number of particles. Ni[t) is a Markov process, which we call 
the "logistic" Markov chain. 

For transition probabilities, we have 

{nbdt + o{dt^) ifj = n + l 

nfidt + -fn'^ / L dt + o{dt'^) ifj = n-l 
o{dt^) otherwise 

If Niit) is large, therefore, there is a left drift, whereas if Ni{t) is small, there 
will be a drift to the right. An important point is the equilibrium point, n^, 
where the rates are equal, that is, 

*2 

bnl = finl + 

which means that 

* 



L{b-fJ.) 
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We will show that as L becomes large, the Markov chain N^lt) tends quickly 
to a neighborhood of and afterward fluctuates randomly around ra^- Our 
goal is the analysis of these fluctuations. 

The logistic Markov chain is a particular case of a birth and death random 
walk on Z\ = {0, 1, . . .}, which has been studied extensively (e.g. [5, 12, 17]). 
The generator for a birth and death random walk, X{t), t > is given by 

= a^ip^x - 1) - (a^- + I3x)'ip{x) + I3^^{x + 1) 

for X > and 

£V(0) = /3o^(l)-/3o^(0). 
For the logistic Markov chain X{t) := N^lt), the transition rates are 

2 

7X 

Px = bx,x > 1, Po = ^, ax = fix -\ — —,x>l 
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We will analyze the logistic Markov chain asymptotically as L — )■ oo. In fact, 
it is more convenient to study a modified logistic chain with transition rates: 



= 6(x + l),x>0, 
Here, the equilibrium point is 



x> 1. 



(3.1) 



(6-/i) + v/(6-/x)2 + 47/L 
27/L 

is equal to the old equilibrium point for large enough L. For convenience, 
when appropriate, we assume that L is such that {h — ij)L/'-^ is an integer. 
Thus, 

_{h- fi)L 

''L — 
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There are two reasons for this change. First, unlike the unmodified chain, 
the modified logistic chain has no absorbing state at x = 0. Second, calcu- 
lations for the modified logistic chain are simpler. Let us stress that as L 
becomes large, the asymptotics are equivalent for both models. 



4 Detailed Analysis of the Logistic Markov 
Chain 

4.1 Hypergeometric Functions 

In this section, we will study the asymptotic properties of the degenerated 
hypergeometric function. For a particular choice of parameters, such func- 
tions will play an important role in the analysis of the logistic Markov chain 
with transition rates (3.1). 

The classical degenerated (or confiuent) hypergeometric function depends 
on two parameters, a and 7, and is given by the Taylor expansion 



az a{a + l)z^ aia + \) . . . (a + n — 1) z'^ 

(h((y ^ r^ = l-| 1 1^ L U -I 1^ 1^ L U 

^ ' ^ 7 1! 7(7 + 1) 2! 7(7 + 1) ... (7 + n - 1) n! 

(see [10], 9.21). $ is an entire function of order 1. We focus on the special 
case when a = 1 and 7 = A, A S> 0, and call such a function 



6 



7{A, z) = 1 + 4 + . . . + -TT-, 7T-^T1 TT + • • • 

^ ' A A{A+ l)...{A + n-l) 

In our case, the general integral representation of $(q;,7,2) ([10], 9.211,1) 
leads to the formula 



1 



-1 

The substitution t = 1 — s leads to 

r-2 



F{A, z) = 2"-^ {A -I) I s^-2e-(i-^)/2 

Jo 

= 2'-^ (A - 1) f s^-\-''l'' ds = ~ r t^'^e-' dt (4.1) 
Jo ^ Jo 

where, in the last step, we used the substitution sz/2 = t. The integral factor 
in (4.1) is the incomplete F- function: 

7(v4- l,z) = [ t^-^e-'dt 
Jo 



Of course, 



t^-^e-* dt 

If A, 2; ^ 1, one can use the Laplace method to obtain asymptotics for 
TlA, z). Putting t = As, we have 



z/A 



A-2 ■ 



7(A -l,z) = A^-i / e-""^'-— '"^^ ds 
Jo 

The critical point of the "phase function" 

^ ^-2, 
$3 = s : — m s 

^ ' A 

(where $'(so) = 0) is Sq = {A — 2)/ A and the asymptotic behavior of 7(^4 — 
l,z) depends on the relationship between z/A and Sq. 

Routine application of the Laplace method leads to the theorem 
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Theorem 4.1. Consider 

If A, z ^ oo and 

I. z < A, \fA/{A - ^) -> 0, then 

//. z > A and \/A/{z — A) ^ 0, then (using Stirling's formula) 

J^iA, z) ~ ~ e/-^+i (^^) (4.2) 

///. 2; = A + /ia/A, h > is a constant, then 



/F. ^ = A — /iv^, h > is a constant, then 

4.2 General Results for ID Random Walks 

Consider the general random walk on Zj'|_ = {0,1,2,...} in continuous time 
and the transition rates from x, ax to the left and Px to the right, where 
(3x > 0,x > and ax > 0,x > 1 and ao = 0. 

The general birth and death random walk with generator 

Cip{x) = axip{x - 1) - {ax + [3x)ip{x) + (3xip{x + 1) 

is ergodic if and only if the following series converges 

S=l + ^ + ^ + ...+ +■■■■ (4.3) 

ai a\a2 a\ . . . On+i 
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In this case, the invariant distribution is given by 

L ai...ax ' 

Let = min{t : x{t) = y}. From the ergodicity of X{t), it follows 
that for any x,y & Z^,x 7^ y, E^Ty is finite. We will use the notations 
ExTy = ETx^y and E^Ty = u{x, y). For fixed y and x > y, u{x, y) satisfies 

Cu{x, ■) = —1, X > y 
uiy,-) = 

One can understand u{x, ■) as lim UL{x,y),x G [y,L], where 

Cul{x, ■) = -1, 
uiiy, ■) = ul{L, ■) = 

Note that for x > y, 

'^x—>-y ~ T~x-^x—l ~l" Tx— 1->X"2 Ty+l— !.y 

Lemma 4.2. ylssnme i/iai t/ie random walk is ergodic. Then 

1 



ETy+l^y 



The second factor has the following interpretation. Put 



(4.5) 



(i.e., 5 = 5o = 1 + ^ + ^ + • • •)• Then, 

7r(x) = lim p(t, -, x) = S~^— ^ . (4.6) 

Sy^ = TTy^'^ is, therefore, the invariant probability for the random walk X{t) 
on [y, 00) with reflection rate f3y at y. 

Considering the cycles between successive transitions y + 1 ^ y and 
y — 7- ?/ + 1 and applying the Law of Large Numbers we deduce that 



vr, 



^/Py 



y 
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and so 

Sy = l+ (3yETy + 

It also follows that ^ 

ET|/+1->J/ = Sy+l. 

CHy+l 

Proof. Let y be fixed. Set 

u{x - y) := Er^^y, x>y. 

Then u{0) = and 

a^u{x - 1) - (a^ + + f3xu{x + 1) = -1. (4.7) 

Set 

A(x) = — u{x — l),x > 1. 

Then 

A(x + 1) = ^A(x) - ^ (4.8) 

Consider a random walk XAr(t) with the same transition rates but on the 
finite interval [0, A^]. We are interested in the first entrance to N. Let 

ro,Ar = min{t : XN{t) = or A^}, Um{x) := ^^{tq^isi)-, and 

AAr(x) := Un{x) — Un{x — 1). 
Then, un satisfies (4.7) with boundary conditions 

M^(0) = un{N) = 

and AAr(x) satisfies (4.8) as well as lim ^n{.x) = A{x), pointwise in x. 

N-^oo 

For the finite chain, 

A^(l) + A;v(2) + . . . + A;v(iV) = (4.9) 
From (4.8) and (4.9), one obtains 

aia2...ak-i. 02 • • • "fc-i 03 . . . Ofc-i 1 



fc-1 
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Note that 

A;v(l) = M7v(l) -Mjv(O) = M^(1). 
We now multiply through by f^f^ to obtain equations for Aiv(l) 

An{1) = AtvIk) H \ h . . . H 

aia2---ttfc-i aia2...ak-i aia2 • • • ai 

This is true for all k, including k = N, and so 

A(l):= hm A^(l) = lfl + ^ + ^ + .. 



Therefore, ETy+i^y is given by 



ETy^.^y = A(y + 1) = ^ ( 1 + + + ...). □ 

Oiy+l \ Oly+2 Oiy+20iy+3 / 

For the logistic chain with 

/3^ = b{x + l) , a^; = X + ^ j , 
we represent Sy, y > in terms of the hypergeometric function F{A,z). 



Oiyj^l Oty+lOiy+2 

1+ * 



1 + 



^+ i (,/^ + — (,/^ + — f— j (4.10) 

, (6L/7)2 



f + y + l (^ + t/ + l)(^ + y + 2 
= F(/iL/7 + y,6L/7). 

This formula, asymptotics from Theorem 4.1, and the obvious relation 

Er^.^y = Er^^^_i + Er^._i^^._2 H h Exy+i^^, for x > y 

will be the basis for the analysis of the first passage times. 
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Let us note also that there is a recurrence formula connecting Er^+i^y 
and Erj/+2->y+i- From (4.5), we obtain 



^'y+l-^y — 1 

Oiy+l CHy+l 



Oiy+2 V (^y+3 OCy^^ay^i 



4.3 Limit theorems for the invariant distribution of the 
logistic Markov chain 

We apply the general results on the ID ergodic random walk on to the 
particular case of the logistic Markov chain, which we defined above as a 
means of studying the mean field Bolker-Pacala process. For the modified 
chain, by the results above, we can obtain a local Central Limit Theorem. 

Theorem 4.3 (Local CLT). Let b > fi. If k = 0{L'^/^), then, for the 
invariant distribution ttl, 

TiL{nl + k) 7==^, 

where a\ = Lb/'-f. 
Proof. From (4.3) we have 
^ ^ b b^ 

S = 1 + 7^ + 7 TTTT ^^TT + - • -4 



/i + 7/L (/i + 7/L)(^ + 27/L) {fi + ^/ L) . . . {fi + n^/ L) 

(4.11) 

And for n G Z^, from (4.4) 

To analyze this series, we consider the position, which we label n*^, such that 
the ratio of the terms ^'^^^^-^-j ~ 1- (If there are two such positions, we take 
to be the larger.) Taking this ratio, we find that 

- 1 < „i < (ii^, 

7 7 
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We obtain the formula 



where A = jj, + = 6 + 0(1/L). Now, we consider 
k ^ 

We integrate the series ln(l+x) = x — + ■ , and take k = O (L^/'^). 



1 

" 2AL~ q' "' 

Thus, 

A: 



,6 + 0(l/L), 

It remains to calculate 7rL(n2), for which we use (4.12). 

We apply the asymptotic formulas for J^{A,z) ((4.2)) to (4.11) to calcu- 
late S: 



^ ^ (1 + ^^Lh) ^ (1 + ^^Lh){2 + /iL/7) ^ ' ' ' 
= ^(/iL/7,6L/7) 

For the product in the denominator, we use the Euler-Maclaurin formula (see 
[20]). First, 
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lini := (/i + l/L) ■ ■ -(/i + n^/L) = (1 + . . . (1 + 



/i^iexp ^ln(l + kuj) 



k=0 



where a; := -y. The Euler-Maclaurin formula gives 



^ In (1 + cu/c) = - / Inxdx + -^^"^{^ + ru) + 0(cu) 

fc=0 ^ 

= (- + r + ^) ln(l + ru)-r + 0{uj). 
00 2 



Thus 



~ ytt'"^ exp 

N Lfe/7+1/2 



1 1 , , 

h r H — ) m(l + rcu) — r 

uj 2 



gL(6-/x)/7 

And so, finally, 



Our large deviations result is 

Theorem 4.4 (Large Deviations). For constant 6 > 0, 
where fiz) := —7 — —-z'^ = / ln(l + x)dx. 
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Proof. Set := n}^{l + S) and k := Applying the formula for gen- 

eral Markov chains (4.4), the calculated invariant probability vr(n2), and the 
Euler-Maclaurin formula (with u := 

7r{N) = Tr{nl) 



ani+i ■ ■ ■ ani+k ^ /27r^L + (^I + " " " (/^ + (^l + 




ln(l + A;cj) - k] 



(4.15) 



Using the expansion ln(l + x) = x — \x'^ + — + . . ., substituting 
in for cu, /c, and n^, and multiplying through, we obtain for the exponent: 

(l + A. + l)ln(l + M-A. = ^-f 4^f^^V. (4.16) 



□ 



4.4 Global Limit Theorems 

A functional Law of Large Numbers and Central Limit Theorem for the 
logistic Markov chain follow directly from theorems of Kurtz [17, 18]. Thus, 
we state them here without proof. 

Define a new stochastic process for the population density, Zi{t) := ^^^^^ . 

Set = !i = 

L 7 

We define the transition function, fii^ii) '■= jPiN^, + j). Thus, 

f ^ = bz J = 1 

fL{z,j)=l t^^^i^ = f^z + jz^ j = -l 

(not needed) j = 

Note that f^^Zjj) does not, in fact, depend on L and we write simply f{z,j). 
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Theorem 4.5 (Functional LLN). As L ^ oo, Zi{t) Z(t) uniformly in 
probability, where Z(t) is a deterministic process, the solution of 

^ = F{Z{t)), Z{0) = z,. (4.17) 

where 

Fi^) ■= ^jfizj) = bz- jiz- -iz^ = -fz{z* - z). 

3 

Equation (4.17) is in fact that of the stochastic logistic model, as discussed 
by Pollett [19]. It has the solution 

Z{t, z) = ^, f\ t > 0. 

z + [z* — z)e ' 

Next, define Gl{z) := ''^^j^fL{z,j) = (6 + + •yz'^. This too does not 

j 

depend on L and we simply write G{z). 

Theorem 4.6 (Functional CLT). // VI(Zl(0) - z*) = (o, the processes 

Cdt) := VliZLit) - Z{t)) 

converge weakly in the space of cadlag functions on any finite time interval 
[0,T] to a Gaussian diffusion ({t) with: 

1) initial value C(0) = Co; 

2) mean 

fF'(Z(u))du 

EC{s) = CoLs := Coeo 

3) variance 

s 

Var(C(s)) = LlJ LZ^G{Z{u))du. 



Suppose, moreover, that F{zq) = 0, i.e., Zq = z* , the equilibrium point. 
Then, Z{t) = zq and({t) is an Ornstein- Uhlenbeck process (OUP) with initial 
value Co J infinitesimal drift 

q := F'{zo) =^i-b 
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and infinitesimal variance 

2b 

a := G{zo) = —{h - ji). 

7 

Thus, C,{t) is normally distributed with mean 
and variance 

Finally, as L — )■ oo, the first exit time tla from n* to n* ± Ay/L will be 
approximately the first exit time ta for the OUP starting at from iA. The 
distribution of ta has been shown by Breiman [3] to be 



P{r^ >t} = ae-^-^^^)* + O (e-2(K^)+'5)t) ^ (4. is) 



such that 



(i) lim u{A) = 

A— >-oo 



(a) lim z/fA) = 00 



[Hi) if is the smallest positive root of 
^ (-2v4^)^ ml 
^ (2A;)! (m-fc)!' 
then = m. 



4.5 First passage time 

As L — 00, 



fc=l fc=l " A;=l ^ 
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Let m = (1 - 5i)n*, Si > 0. Then 

n* n* ^ 

A;=rii + 1 fc=ni+l " (4.19) 

'S'ni+1 A ni + 1 ni + 1 2 , \ 1 



1 TTtPi T^Pi + • • • 



Ofni+i V ^1 + 2 ni + 3 / 1 - pi 

where pi := 1 — (1 — f 

Next, we analyze the first passage time u{x) to {ni,n2} for x G [ni,n2], 
where rii = (1 — (5i)n*, n2 = (1 + ^2)^^*5 5i, ^2 > 0. u must satisfy 

£m = -1, (4.20) 

as well as the boundary conditions u{n\) = u{n2) = 0. A particular solution 
to (4.20) is ipii^) '■= Er^_>.o- Then, 

u = 'lpi + Ci + C21IJ2, 

where ci and C2 are constants and ip2 satisfies Cip2 = 0. 

We choose ip2 to satisfy 4'2{n*) = and 4'2{n* + !) = !. This gives 

l3n*Pn*-l I _j_ ^n*---fe + l \ if X < 71* 



M^)= { ' iix = (4.21) 

1 + + . ■ . + ifx>n* 

Pn* + 1 Pn*+l'"Pa:-l 

Let us find the asymptotics of tp2{x) for a; = ni = n*{l — 6i), x = 77.2 = 
77* (1 + ^2)- We will calculate up to a constant factor and use the standard 
notation, a„ x 6„ to mean 

0<Ci<--^<C2<OO 
On 

for appropriate constants Ci and C2. 
The last term 

Pn* + 1 ■ ■ ■ Pn2-1 

gives the main contribution to 4^2(^2) ■ The terms An2-k asymptotically form 
a geometric progression with the ratio 

1 

~ l + 62{l-fi/b) ^ ^ 
18 



(because 



«n2-fc ^ /i(r22 -k) + l{n2 - kf ^ ^ 7 

/3n2-fc ^('^2 - /C + 1) 

= i + 5,ri-^) + ori).) 



6 + 6L^^ + ^^^^ 



This means that, for a; := 



z/'2('^2) ~ ^n.2-l(l + P2 + P2 + • • • 



1 - P2 



Moreover, using the Euler-Maclaurin formula (compare (4.13)) 



/ 52(1-^) 

1 



hri* /52n* \ (if 

An2-i X ]^(l+/ccj) = exp ^ln(l + ku) j x exp — / 

k=0 \k=0 J Q 



ln(l + x)dx 



J 



and so 



■02(^2) X exp 



/ 52(1-^) 

1 



\ 



\ 



ln(l + x)dx 



I 



Similar calculations for '^i(n\) yield 



^2(^1) X - exp 



V 



-— / ln(l — x)dx 
^ J 



J 



It is convenient to introduce some sort of "symmetry" of the logistic 



Markov chain with respect to the equilibrium point n* 
simplest way to do so is to assume that 

Pn*{^L(^[ni,n2]) = ^^2} ~ ^• 

This means that ip2{n2) ^ —'ip2i'ni), i.e. 

52(1-^) 



(b-li)L 



The 



ln(l — x)dx 



ln(l + x)dx. 
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The last equation uniquely determines 62 as a function of < (5i < 1. 

With this symmetry, we can determine Er„._i.{„j^„2}. Consider again the 
problem 

Cu = -1, u{ni) = u{n2) = 0. (4.22) 

We modify the previous particular solution and define ipi^x) := Er^^^m for 
X > rii. Then, 

u{x) := ijji{x) + Ci + C2ip2{x) 

satisfies 4.22 for some constants ci and C2- Using ip2{n2) x —ip2{ni) and 
'^i('^i) = 0, we obtain ci = — |'?/^i(n2). From 4.19 and using I of Theorem 4.1 
we have that 

V;i(n2) = Er„2^„, ^J!±±}: 



Q!ni + l(l - Pi) 

A/SvFyexp ^^Llnpi + 5i{l - lnpi)n"' 
(6-p)(l-(5i)(l-pi)v/5;^ " 
exp (^^Llnpi + (5i(l -lnpi)n*^ 

7l ^ 



And so, finally 



^ exp(^Llnpi + 5i(l-lnpi)n*) ^^'^S) 

^ 7l • 

This result tells us that arrival at is extremely improbable and, there- 
fore, we can essentially ignore it. For the central limit theorem and large 
deviation results of the previous sections, therefore, we do not need to as- 
sume ergodicity, e.g., by adjusting the /3 parameter to prevent absorption 
at 0. This also suggests that the expected recurrence time to a point k, 
Tfc ~ 7r{k)~^, with 7r{k) the invariant probability of state k. In particular, the 
expected recurrence time to the equilibrium point will be 0{\/L) and to 
a point nl + SL will be yfLe^^^\ 
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